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Abstract. This is an introduction to elementary decoherence theory as it is typically 
applied to superconducting qubits. 

Abbreviations: SQUID - superconducting quantum interference device; qubit - quan- 
tum bit; TSS - two state system 
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1. Introduction 

. The transition from quantum to classical physics, now known as decoher- 
\Q \ ence, has intrigued physicists since the formulation of quantum mechanics 
O . (Giulini et al., 1996; Leggett, 2002; Peres, 1993; Feynman and Vernon, 
1963; Zurek, 1993). It has been put into the poignant Schrodinger cat para- 
dox (Schrodinger, 1935) and was considered an open fundamental question 
for a long time. 
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In this chapter, we study the theory of decoherence as it is apphed 
to superconducting qubits. The foundations of the methodology used are 
rather general results of quantum statistical physics and resemble those 
applied to chemical physics, nuclear magnetic resonance, optics, and other 
condensed matter systems (Weiss, 1999). All these realizations introduce 
their subtleties — typical couplings, temperatures, properties of the cor- 
relation functions. We will in the following largely stick to effective spin 
notation in order to emphasize this universality, still taking most of the 
examples from superconducting decoherence. This paper is based on lec- 
tures 2 and 3 of the NATO-ASI on "Manipulating quantum coherence in 
superconductors and semiconductors" in Cluj-Napoca, Romania, 2005. It 
is not intended to be a review summarizing the main papers in the field. 
Rather, it is an (almost) self-contained introduction to some of the relevant 
techniques, aimed to be accessible to researchers and graduate students 
with a knowledge of quantum mechanics ( Cohen- Tannoudji et al., 1992) 
and statistical physics (Landau and Lifshitz, 1984) on the level of a first 
graduate course. So much of the material here is not new and most certainly 
known to more experienced researchers, however, we felt a lack of a single 
reference which allows newcomers to get started without excessive overhead. 
References have largely been chosen for the aid they provide in learning and 
teaching the subject, rather than importance and achievement. 

1.1. BASIC NOTIONS OF DECOHERENCE 

The mechanisms of decoherence are usually related to those of energy dis- 
sipation. In particular, decoherence is irreversible. If we take as an example 
a pure superposition state 



and compare it to the corresponding classical mixture leading to the same 
expectation value of ct^ 



we can see that the von-Neumann entropy p = —ksTi [plog p] rises from 
•Spure = to Smix = 1^2. Hencc, dccoherencc taking Ppurc to pmix creates 
entropy and is irreversible. 

Quantum mechanics, on the other hand, is always reversible. It can 
be shown, that any isolated quantum system is described by the Liouville 
von-Neumann equation 




(1) 




(2) 



ihp = [H, p] 



(3) 
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which conserves entropy. Indeed, also the CPT theorem of relativistic quan- 
tum mechanics (Sakurai, 1967) states, that for each quantum system it is 
possible to find a counterpart (with inversed parity and charge) whose 
time arrow runs backwards. The apparent contradiction between microre- 
versibility — reversibility of the laws of quantum physics described by 
Schrodinger's equation — and macro-irreversibility is a problem at the 
foundation of statistical thermodynamics. We also remark that the La- 
grangian formalism (Landau and Lifshitz, 1982) which was used as the 
starting point in the previous chapter of this book (Geller et al., 2006) does 
not even accomodate friction on a classical level without artificial and in 
general non-quantizable additions. 

1.1.1. Heat baths and quantum Brownian motion 

The standard way out of this dilemma is to introduce a continuum of ad- 
ditional degrees of freedom acting as a heat bath for the quantum system 
under consideration (Feynman and Vernon, 1963; Caldeira and Leggett, 
1981; Caldeira and Leggett, 1983). The complete system is fully quantum- 
coherent and can be described by equation 3. However, the heat bath 
contains unobserved degrees of freedom which have to be integrated out 
to obtain the reduced system; the reduced system is the original quantum 
system which does not contain the bath explicitly, but whose dynamics are 
influenced by the bath. The dynamics of the reduced system now show both 
dissipation (energy exchange with the heat bath) and decoherence (loss of 
quantum information to the heat bath). Another view on this is that any 
finite combined quantum system shows dynamics which are periodic in 
time. The typical periods are given by the inverse level splittings of the 
system. Thus, a continuous heat bath shows periodicity and reversibility 
only on an infinite, physically unobservable time scale. 

A standard example, taken from Ref. (Ingold, 1998), of irreversibility 
in both classical and quantum mechanics is (quantum) Brownian motion 
(QBM), which we will now describe in the one-dimensional case. The under- 
lying Hamiltonian of a single particle in an oscillator bath has the general 
structure 

H = Hs + H,b + Hb + H,. (4) 

Here, the system Hamiltonian Hs describes an undamped particle of mass 
M in a scalar potential, Hg = ^ + V{q). Hi, describes a bath of harmonic 

oscillators. Hi = J2i + ■ The coupling between these two 

components is bilinear, Hsb = —qJ2i(^i^i- If this were all, the effective 
potential seen by the particle would be altered even on the classical level, 
as will become more obvious later on. Thus, we have to add a counter term 
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which does not act on the bath, He = q J2i o — Adding this counterterm 
gives the Hamiltonian the fohowing intuitive form 

indicating that the bath oscillators can be viewed as attached to the particle 
by springs. Here, we have introduced sets of new parameters, q, coi, and 
rui which need to be adjusted to the system of interest. This aspect will be 
discussed later on. We treat this system now using the Heisenberg equation 
of motion 

ihd{t) = [0{t),H] (6) 

for the operators q, P, Xj, and pi, which (as a mathematical consequence 
of the correspondence principle) coincide with the classical equations of 
motion. The bath oscillators see the qubit acting as an external force 

Xi+ujfxi = —q{t). (7) 
rrii 

This equation of motion can be solved by variation of constants, which can 
be found in textbooks on differential equations such as (Zill, 2000) 

Xi{t) = Xi{0) cos uJit -\ — sinwjtH ^ / dt' s\nuji{t — t')q{t') (8) 

rriiUOi ruiUof Jq 

Analogously, we find the equation of motion for the particle 

^• = -?^-E--^-9E^- (9) 

oq ^ rrii ^ rriiUJf 

Substituting eq. 8 into eq. 9 eliminates the bath coordinates up to the initial 
condition 

Mq = -^-V-^ / dt' sin (jji(t - t')q(t') 
dq ^i^i -^0 

+ y Q fxi(0)cosa;,t + ^^sinwit) - (10) 
^ V rriiUJi J Y rriiLof 

We now integrate by parts and get a convolution of the velocity plus bound- 
ary terms, one of which shifts the origin of the initial position, the other 
cancels the counterterm (indicating, that without the counterterm we would 
obtain a potential renormalization) . The result has the compact form 

dV /■* 

+ — + dt'^it - t')q{t') = m- (11) 
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This structure is identified as a Langevin equation with memory friction. If 
interpreted classically, this is the equation of motion of a Brownian particle 
- a light particle in a fluctuating medium. In the quantum limit, we have 
to read g, Xj and the derived quantity ^ as operators. We see both sides of 
open system dynamics — Dissipation encoded in the damping kernel 7 and 
decoherence encoded in the noise term ^. We can express 7 as 

jU) = / ^— ^ COS U!it= / J(u))cosu!t (12) 

^ rriiLof Jo to 

where we have introduced the spectral density of bath modes 

j[^) = y -i^i) (13) 

which is the only quantifier necessary to describe the information encoded 
in the distribution of the rrii, uji, and Cj. The right hand side of eq. 11 is a 
noise term and reads 



e(t) = E' 



Xii^j 9<?(Uj COS uJit-\ smwjt 

miujf I niiUJi 



(14) 



This crucially depends on the initial condition of the bath. If we assume 
that the bath is initially equilibrated around the initial position q{{)) of 
the particle, we can show, using the standard quantum-statistics of the 
simple harmonic oscillator, that the noise is unbiased, (^(t)) = 0, and its 
correlation function is given by 

K{t) = (^(t)C(O)) = / (hjJ{Lo) [cos ujt {2n{hLo) + sin uot] (15) 



where n is the Bose function, n{huj) = (g^"^/*^^ — 1)~^, and 2n{fujj) + 1 = 
coth 2kBT ) • Here and henceforth, angular brackets around an operator 
indicate the quantum-statistical average, (O) = Tr(pO) with p being the 
appropriate density matrix. We will get back to the topic of the initial 
condition in section 3.1.2 of this chapter. 

The noise described by (t) is the quantum noise of the bath. In partic- 
ular, the correlation function is time-translation invariant, 

K{t) = {at + T)ar)) (16) 

but not symmetric 



K{-t) = mm) = K\t) / K{t). 



(17) 
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which reflects the fact that ^ as defined in eq. 14 is a time-dependent 
operator which does generally not commute at two different times. Ex- 
phcitly, the imaginary part of K(t) changes its sign under time reversal. 
Indeed, if the derivation of eq. 15 is done explicitly, one directly sees that 
it originates from the finite commutator. Moreover, we can observe that 
at T ^ w we have 2n + 1 — > 2kBT/huj ^ 1, thus the integral in eq, 15 is 
dominated by the symmetric real part now describing purely thermal noise. 
At any temperature, the symmetrized semiclassical spectral noise power in 
frequency space reads 



where (. . .)^ means averaging and Fourier transforming. This quantity con- 
tains a sign of the quantum nature of noise. Unlike classical noise, it does 
not disappear at low temperatures T ^ Tiuj/kB, but saturates to a finite 
value set by the zero-point fluctuations, whereas at high temperature we 
recover thermal noise. Note, that the same crossover temperature dictates 
the asymmetry in eq. (17). Both observations together can be identified with 
the fact, that zero-point fluctuations only allow for emission of energy, not 
absorption, as will be detailed in a later section of this chapter. 

Our approach in this chapter is phenomenological. The main param- 
eter of our model is the spectral density J{lo). We will show in sections 
1.1.3, 2.1.1 and 2.1.2 how J{io) can be derived explicitly for Josephson 
junction circuits. Oscillator baths accurately model numerous other situa- 
tions. Decoherence induced by phonons in quantum dot systems allows to 
directly identify the phonons as the bath oscillators (Brandes and Kramer, 
1999; Storcz et al., 2005a), whereas in the case of electric noise from resistors 
or cotunneling in dots (Hartmann and Wilhelm, 2004) it is less obvious 
— the Bosons are electron-hole excitations, which turn out to have the 
commutation relation of hard-core bosons (von Delft and Schoeller, 1998) 
with the hard-core term being of little effect in the limits of interest (Weiss, 



Going back to our phenomenology, we introduce the most important 
case of an Ohmic bath 



Here, 7 is a constant of dimension frequency and / is a high-frequency cutoff 
function providing f[x) ~ 1 at a:; < 1 and / ^ at x > 1. Popular choices 
include the hard cutoff, f{x) = 6{1 — x), exponential cutoff, f{x) = e~^, 
and the Drude cutoff f{x) = j^^- We will see in section 1.1.3 that the 
Drude cutoff plays a significant role in finite electrical circuits, so we chose 
it here for illustration purposes. In this case, the damping kernel reduces 



(18) 



1999). 




(19) 
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to 



For uJc 




damping with damping constant 7, 7(t) = 7(^(t). Here, "classical damping" 
alludes to the damping of particle motion in fluid or of charge transport 
in a resistor (thus the name Ohmic, see also section 1.1.3). With finite 
cjc, the Ohmic models leads to classical, linear friction proportional to the 
velocity, smeared out over a memory time set by the inverse cutoff frequency 
defining a correlation time tc = ^c^- On the other hand, as it turns out in 
the analysis of the model e.g. in section 2.2.2, an infinite cutoff always leads 
to unphysical divergencies. Examples will be given later on. All examples 
from the class of superconducting qubits have a natural ultraviolet cutoff 
set by an appropriate 1/RC or R/L with R, L, and C being characteristic 
resistances, inductances, and capacitances of the circuit, respectively. Note, 
that parts of the open quantum systems literature do not make this last 
observation. 

We will not dwell on methods of solution of the quantum Langevin equa- 
tion, as the focus of this work is the decoherence of qubit systems. Methods 
include the associated Fokker-Planck equation, path integrals, and quantum 
trajectory simulations. The quantum Langevin equation finds application in 
the theory of quantum decay in chemical reactions, the dissipative harmonic 
oscillator, and the decoherence of double-slit experiments. 

1.1.2. How general are oscillator baths? 

Even though the model introduced looks quite artificial and specific, it 
applies to a broad range of systems. The model essentially applies as long 
as the heat bath can be treated within linear response theory, meaning that 
it is essentially infinite (i.e. cannot be exhausted), has a regular spectrum, 
and is in thermal equilibrium. We outline the requirement of only weakly 
perturbing the system, i.e. of linear response theory (Kubo et al., 1991). 
The derivation is rather sketchy and just states the main results because 
this methodology will not be directly used later on. Introductions can be 



found e.g. in Ref. (Kubo et al., 1991; Ingold, 1998; Callen and Welton, 



In linear response theory, we start from a Hamiltonian Hq of the oscil- 
lator bath which is perturbed by an external force F coupling to a bath 
operator Q, 



where the perturbation must be weak enough to be treated to lowest order. 
It is a result of linear response theory that the system responds by a shift 



1951) 



H = Ho-FQ 
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of Q (taken in the Heisenberg picture) according to 

{m))^ = x{^)F{u) (22) 

where the susceptibiUty x can be computed to lowest order as the correla- 
tion function 

X = {Qmimt))^ (23) 

computed in thermal equilibrium. We can split the correlation function 
into real and imaginary parts X = x' + ^x"- The real part determines the 
fluctuations, i.e. 

1 {5Q(t)5Q{d) + 5Q{Q)5Q{t))^ = x' (24) 

whereas the imaginary part determines the energy dissipation 

{E{t))^=u:x"\F\\ (25) 

together with the equations 12 and 18 tracing both damping and noise 
back to a single function x constitute the famous fluctuation-dissipation 
theorem (Callen and Welton, 1951), a generalization of the Einstein relation 
in diffusion. 

In this very successful approach we have characterized the distribution 
of the observable Q close to thermal equilibrium by its two-point correlation 
function alone. This is a manifestation of the fact that its distribution, 
following the central limit theorem is Gaussian, i.e. can be characterized 
by two numbers only: mean and standard deviation. Oscillator baths pro- 
vide exactly these ingredients: by properly chosing J{oo) they can be fully 
adjusted to any xi}^)^ aiid all higher correlation functions — correlation 
functions involving more than two operators — can also be expressed 
through J hence do not contain any independent piece of information. 

This underpins the initial statement that oscillator baths can describe a 
broad range of environments, including those composed of Fermions and not 
Bosons, such as a resistor. As explained in section 1.1.1, the oscillators are 
introduced artifically — on purely statistical grounds as a tool to describe 
fluctuations and response — and can only sometimes be directly identified 
with a physical entity. 

There are still a number of environments where the mapping on an 
oscillator bath is in general not correct. These include i) baths of uncoupled 
spins (e.g. nuclear spins), which are not too big and can easily saturate, 
i.e. explore the full finite capacity of their bounded energy spectrum ii) 
shot noise, which is not in thermal equilibrium iii) nonlinear electrical 
circuits such as many Josephson circuits and iv) in most cases 1// noise. 
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whose microscopic explanation either hints at non-Gaussian (spin-hke) or 
nonequihbrium sources as discussed in section 3.1.3. 

1.1.3. Oscillator bath models for Josephson junction devices 
We have now learned two approaches to characterize the oscillator bath: 
through noise, and through friction. We will now apply the characterization 
by friction to a simple Josephson circuit with Josephson energy Ej, junction 
capacitance Cj and arbitrary shunt admittance in parallel, all biased by an 
external current Ib- We are extending the method presented in the previous 
chapter (Geller et al., 2006) to include the admittance. We start with the 
elementary case of a constant conductance, Y{lu) = G. The total current 
splits up into the three elements as 

lB = IcSi^ct> + C^'<j) + G^^. (26) 
zvr zvr 

Reordering terms, we can cast this into the shape of Newton's equation for 
a particle with coordinate (j). 

Here, we have multiplied the equation by another to ensure proper 

dimensions of the potential energy 

V{(t)) = -lBt^(t> + Ej{l-cos(t>) (28) 

ZTT 

where we have introduced the Josephson energy Ej = Io^o/^tt. This ex- 
pression can be readily compared to eq. 11. We see that the friction term 
has no memory, i.e. ^{t) oc 5{t), and using the results of section 1.1.1 we 
can infer that J(w) = G{^o/2tt)u>, i.e. an Ohmic resistor leads naturally 
to an Ohmic spectral density as mentioned before. Note that this has no 
cutoff, but any model of an Ohmic resistor leads to reactive behavior at 
high frequencies. 

We see that we missed the noise term on the right, which would represent 
current noise originating in G and which would have to be included in a 
more sophisticated circuit analysis which careful engineers would do. By 
applying the fluctuation dissipation theorem to 7 we can add on the proper 
noise term, whose correlation function is given by equation (15) — or we 
can simply use this equation with the J{uj) obtained. 

We want to generalize this system now to an arbitrary shunt admittance 
Y(ij). For that, it comes in handy to work in Fourier space and we denote 
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the Fourier transform by T . Analogous to eq. (27), we can find the following 
expression 

-.^c(^)%....m(|)%..(|^)^o. 

We have to remember that the damping Kernel 7 is the Fourier cosine 
transform of J(u;)/a;, which also implies that it is a real valued function. 
We can split Y into real (dissipative) and imaginary (reactive) parts Y = 
Yd + iYr- For any finite electrical circuit, Y^ is always an even and Yr always 
an odd function of frequency. All this allows us to rewrite eq. 29 

-.^c(|)%-.KH(£)%...w(|)%..(-).„. 

(30) 

Thus, the general expression for the spectral density reads J{uj) = ujYd = 
a;Rey(ti;), i.e. it is controlled by the dissipative component of Y{ijj) alone. 
There is a new term containing the reactive component Y^ which modifies 
the non-dissipative part of the dynamics and can lead e.g. to mass or 
potential renormalization, or something more complicated. Comparing this 
result to the structure of the susceptibility % in the discussion of section 
1.1.2 it looks like the real and imaginary part have changed their role and 
there is an extra factor of w. This is due to the fact that Y links I and V ^ 
whereas the energy-valued perturbation term in the sense of section 1.1.2 
is QV . This aspects adds a time-derivative Y = x which leads to a factor 
iuj in Fourier space. 

This last result can be illustrated by a few examples. If Y[u)) = G<I>o/27r, 
we recover the previous equation (27). If the shunt is a capacitor C^, we 
have Y{ijj) = iioCs and we get from eq. (30) the equation of motion of 
a particle with larger mass, parameterized by a total capacitance Ctot = 
Cj + Cg. On the other hand, if the shunt is an inductance Lg? we obtain 
Y{ijj) = {iujLs)~^, leading to a new contribution to the potential originating 
from the inductive energy 

and no damping term. Finally, let us consider the elementary mixed case 
of a shunt consisting of a resistor Rg and a capacitor Cg in series. We 
find Y{io) = ^^ZrIcs ^^i^^ 

can be broken into a damping part which is 
supressed below a rolloff frequency uJr = {RC)~^, 1^-^ ^^^^ ^ reac- 

tive part which responds capacitively below that rolloff, Yr = ^"^C* j-j-^^^y^ . 
As the rolloffs are very soft, there is no straightforward mapping onto a 
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very simple model and we have to accept that the dynamics get more 
complicated and contain a frequency-dependent mass and friction as well 
as time-correlated noise, all of which gives rise to rich physics (Robertson 
et al., 2005). 

2. Single qubit decoherence 

2.1. TWO-STATE OSCILLATOR BATH MODELS 

In the previous section, we introduced the notion of an oscillator bath envi- 
ronment for continuous systems including biased Josephson junctions. We 
derived quantum Langevin equation demonstrating the analogy to classical 
dissipative motion, but did not describe how to solve them. In fact, solving 
these equations in all generality is extremely hard in the quantum limit, 
thus a restriction of generality is sought. For our two-state systems (TSS) of 
interest, qubits, we are specifically interested in the case where the potential 
in the Hamiltonian of eq. 5 forms a double well with exactly one bound 
state per minimum, tunnel-coupled to each other and well separated from 
the higher excited levels, (Geller et al., 2006). When we also concentrate on 
the low-energy dynamics, we can replace the particle coordinate q by qqcTz 
and the Hamiltonian reads 



H = i:'yz + ^<yx + ^Y. + "1) + E + 1/2), (32) 



where e is the energy bias and A is the tunnel splitting. This is the fa- 
mous Spin-Boson Hamiltonian (Leggett et al., 1987; Weiss, 1999). We have 
dropped the counterterm, which is oc in the continuous limit and, due 
to g = ibgo is constant in the two-state case. The spectral density is 
constructed out of the J{ijo) in the continuous limit 



The Spin-Boson Hamiltonian, eq. (33) is more general than the truncation 
of the energy spectrum in a double-well potential may suggest. In fact, 
it can be derived by an alternative procedure which performs the two- 
state approximation first (or departs from a two-state Hamiltonian without 
asking for its origin) and then characterizes the bath. The oscillator bath 
approximation holds under the same conditions explained in section 1.1.2 
The Spin-Boson model makes the assumption, that each oscillator couples 
to the same observable of the TSS which can always be labelled Uz- This is 



A 




(33) 
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a restrictive assumption which is not necessarily true for aU reahzations of 
a dissipative two-state system. 

As the two-state counterpart to classical friction used in the continuous 
case is not straightforward to determine, the environmental spectrum is 
computed from the semiclassical noise of the environment, following the 
prescription that, if we rewrite eq. 32 in the interaction picture with respect 
to the bath as 

Hi = + ^a., + ^ u,{a\a, + 1/2) (34) 

i 

we can identify 5e for any physical model mapping on the Spin-Boson model 
as 

i {5e{t)5e{Q) + 5e{Q)5e{t))^ = Jtss{co) coth (^^) (35) 
An application of this procedure will be presented in the next subsection. 

2.1.1. Characterization of qubit environments through noise 
A standard application of the characterization of the environment is the 
description of control electronics of relatively modest complexity, attached 
to a flux qubit. We look at the definite example shown in 2.1.1. It shows 
a simplified model of the microwave leads providing the control of a flux 
qubits. The microwaves inductively couple to the sample by a mutual in- 
ductance M between the qubit and a coil with self-inductance L. These 
leads are mounted in the cold part of the cryostat, usually on the qubit 
chip, and are connected to the outside world by a coaxial line which almost 
inevitably has an impedance oi Z = 5017. That impedance provides — in 
light of the discussion in the previous section — a significant source of 
damping and decoherence. As a design element, one can put two resistors 
of size R close to the coil. 

The environmental noise is easily described by the Nyquist noise (Callen 
and Welton, 1951) of the voltage V between the arms of the circuit, see 
figure 2.1.1. The Johnson-Nyquist formula gives the voltage noise 

Sv = l {V{t)V{0) + V{0)V{t)) = hioReZ^s coth (^^) (36) 

where Z^s is the effective impedance between the arms, here of a parallel 
setup of a resistor and an inductor 

and Loff is the total impedance of the coupled set of conductors as seen 
from the circuit. For microwave leads, the total inductance is dominated 
by the self-inductance of the coil, hence L^s ~ L. 
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X XX 



Figure 1. Typical on-chip electromagnetic environment of a superconducting flux qubit, 
consisting of the flux control coil with self-inductance L, mutual inductance M to the 
qubit, shunt impedance Z and on-chip decoupling resistors R. 



We need to convert the voltage noise into energy level noise of the qubit. 
A voltage fluctuation 6V leads to a current fluctuation in the coil following 

81 = 5V/iujL. (38) 

The current noise produces flux noise through the qubit loop 

M 

5$ = MSI = —5V (39) 
iujL 

which converts into energy bias noise following 

MI 

5e = Is6^ = — ^6V (40) 
iojL 

with Is being the circulating current in the potential minima of the qubit. 
Thus, the energy level correlation function reads 

which allows us to express the spectral density through the impedance as 

/MT \ ^ 

J{uj) = hio ReZ^siu). (42) 

V iojL J 

With the specific circuit shown in figure 2.1.1, we find that the environment 
is Ohmic with a Drude cutoff 
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with iOc = L/R and a 



. Thus, we find a simple method to engineer 



h{Z+2R) 



the decoherence properties of thw circuit with our goal being to reduce 
J{lo) by decoupling the device from the shunt Z. The method of choice is 
to put large resistors R on chip. Their size wih ultimately be limited by the 
necessity of cooling them to cryogenic temperatures. The friction method 
introduced earlier, section 1.1.3 , leads to the same result. 

2.1.2. Linearization of nonlinear environments 

In general, nonlinear environments important for qubit devices can also 
be identified. In superconducting devices, these include electronic environ- 
ments which in addition to the linear circuit elements discussed in the 
previous section, also contain Josephson junctions. In general, such envi- 
ronments cannot be described by oscillator bath models, whose response 
would be strictly linear. Here, we want to concentrate on the case of a 
nonlinear environment — a SQUID detector — in the regime of small signal 
response, i.e. in a regime where it can be linearized. This linearization can be 
illustrated by the concept of Josephson inductance. Let us remind ourselves, 
that a linear inductor is defined through the following current-flux relation 



whereas the small flux-signal response of a Josephson junction can be 
approximated as 



where we have split the flux into its average <!> and small deviations (5$ and 
have introduced the Josephson inductance Lj = $o/2vr Ic cos 0. Thus, the 
small-signal response is inductive. 

We would now like to demostrate this idea on the example of a DC- 
SQUID detector inductively coupled to the qubit, see flg. 2.1.2. 

In the first stage, we again need to find the voltage noise between the 
branches of the circuit. This is given by eq. (36) with the appropriate 
inductance calculated from the cicruit shown in the lower panel of fig. 2.1.2, 
= R~^ + iojC + {iujLj)~^ . This is the impedance of an LC resonator 
with damping. The conversion into energy level noise goes along similar 
lines as before, incorporating the SQUID equations as described here and 
in standard literature (Tinkham, 1996; Clarke and Braginski, 2004). 

The DC-SQUID is a parallel setup of two Josephson junctions 1 and 2, 
which for simplicity are assumed to be identical. The total current flowing 
through the device is 




(44) 



/ = 



sin 27r 




(45) 



Ib = ^c(sin (f)i + sin 4)2) = 2Ic cos((5(^/2) sin 



(46) 
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Figure 2. Upper panel: DC-SQUID readout circuit consisting of the actual SQUID, a 
shunt capacitor, and a voltmeter with an unavoidable resistor. Lower panel: Linearized 
circuit used for the noise calculation 



where we have introduced (j) = (01 + (/>2)/2 and Scj) = 4>i — (j)2- Now we need 
to remember that the phases <j)i are connected to the Schrodinger equa- 
tion for the superconducting condensate. Thus, an elementary calculation 
(Tinkham, 1996; Clarke and Braginski, 2004) leads to 

^> 

50 = 27r— mod27r (47) 

where $ is the total magnetic flux through the loop. This is identical to the 
flux applied externally using a biasing coil plus the qubit flux as we neglect 
self-inductance. Thus, for the bias current Ib the DC-SQUID acts like a tun- 
able Josephson junction with a critical current /c,eff(^) = 2/c| cos(7r$/<I>o)|- 
Thus, we can translate voltage fluctuations into phase fluctuations as 



(48) 



16 F.K. WILHELM, M.J. STORCZ, U. HARTMANN, AND M. GELLER 



The qubit is coupling to the magnetic flux which — assuming a symmetric 
SQUID geometry - is coupled only to the circulating current 

/circ = /c(sin0i - sin</)2)/2 = IcCos{<p) sin — . (49) 

We can now express its fluctuations through the fluctuations of (p 

7r<I> - - Ib vr<I> 
<^-^circ = —Ic sin — — sm{(l))S(f) = — tan -^Scj) (50) 

where in the last step we have used eq. 46. With the remaining steps 
analogous to the previous section, we obtain 

/ /r 27r 7r<l>\^ 
J{uj) = fiuj Mls-^—t&n — ReZcff. (51) 

Here, Z^s is the impedance of the linearized circuit shown in the bottom 
panel of fig. 2.1.2. This result reveals a few remarkable features. Most 
prominently, it shows that J{uj) can be tuned by shifting the working point 
of the linearization through changing the bias current Ib- In particular, 
J{oj) can be set to zero by chosing Ib = 0. The origin of this decoupling 
can be seen in eq. 50, which connects the bias current noise to the circulating 
current noise. The physical reason for this is, that in the absence of a bias 
current the setup is fully symmetric — any noise from the external circuitry 
splits into equal amounts on the branches of the loop and thus does not 
lead to flux noise. For a detector, this is a highly desired property. It allows 
to switch the detector off completely. When we do not bias, we have (for 
the traditional switching current measurement) no senitivity and with it no 
backaction. This means, that if the device is really highly symmetric, one 
can push this device to the strong measurement regime while still being able 
to operate in the "off' state of the detector. This effect has been predicted 
in Refs. (van der Wal et al., 2003; Wilhelm et al., 2003). Experimentally, it 
was first observed that the decoupled point was far from zero bias due to 
a fabrication issue (Burkard et al., 2005), which was later solved such that 
our prediction has indeed been verified (Bertet et al., 2005a). 

2.1.3. The Block equation 

So far, we have discussed the characterization of the environment at length. 
We did not specify how to describe the qubit dynamics under its influence. 
For a continuous system, we have derived the quantum Langevin equa- 
tion (11). Even though this eqution looks straightforward, solving it for 
potentials others than the harmonic oscillator is difficult without further 
approximations. We will now show first how to describe decoherence in 



SUPERCONDUCTING QUBITS II: DECOHERENCE 17 



a phenomenological way and then discuss how to reconcile microscopic 
modehing with the Bloch equation. 

For describing the decoherence of a qubit we have to use the density 
matrix formahsm. which can describe pure as well as mixed states. In 
the case of a qubit with a two-dimensional Hilbert space, we can fully 
parameterize the density matrix by its three spin projections Si = Ti{pai), 
i = x,y, z as 

P=l(^l + J2S^^i^ (52) 

where the cij are Pauli matrices. This notation is inspired by spin reso- 
nance and is applicable to any two-state system including those realized 
in superconducting qubits. We can take the analogy further and use the 
typical NMR notation with a strong static magnetic field Bz{t) applied 
in one direction identified as the z-direction and a small AC field, Bx{t) 
and By(t) in the xy-plane. In that case, there is clearly a preferred-axis 
symmetry and two distinct relaxation rates, the longitudinal rate 1/Ti and 
the transversal rate I/T2 can be introduced phenomenologically to yield 

= j{B X S). - ^1-^ (53) 

S,/y = j{B X 5),/, - ^ (54) 

where we have introduced the equilibrium spin projection Sz,eq and the 
spin vector S = {Sx, Sy, S^)'^ . Note that the coherent part of the time 
evolution is still present. It enters the Bloch equation via the Hamiltonian, 
decomposed into Pauli matrices as H = —^B ■ S. This spin notation is also 
useful for superconducting qubits, even though the three components usu- 
ally depend very distinct observables such as charge, flux, and current. This 
parameterization leads to the practical visualization of the state and the 
Hamiltonian as a point and an axis in three-dimensional space respectively. 
The free evolution of the qubit then corresponds to Larmor precession 
around the magnetic field. The pure states of the spin have 5^ = 1 and 
are hence on a unit sphere, the Bloch sphere, whereas the mixed states are 
inside the sphere — in the Bloch ball. 

The rates are also readily interpreted in physical terms. As the large 
static field points in the z-direction in our setting, the energy dissipation is 
given as 

^ = -iBjz (55) 
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and hence its irreversible part is given through l/Ti. On the other hand, 
the purity (or hnearized entropy) P = Trp^ = 1/4 + Sf decays as 

P = 2Y^SA = -^l±^ - - ^-•-^) (56) 

thus all rates contribute to decoherence. Note, that at low temperatures 
Sz,eq — > 1 SO the Ti-term in general augments the purity and reestablishes 
coherence. This can be understood as the system approaches the ground 
state, which is a pure state. In this light, it needs to be imposed that P < 1 
as otherwise the density matrix has negative eigenvalues. This enforces 
T2 < 2Ti. 



2.2. SOLUTIONS OF THE BLOCK EQUATION AND SPECTROSCOPY 

The rates shown in the Bloch equation are also related to typical spectro- 
copic parameters (Abragam, 1983; Goorden and Wilhelm, 2003). We chose 
a rotating driving field 

Bx = {ujRh)cosujt (57) 
By = (wR/7)sinwt. (58) 

In spectroscopy, we are asking for the steady state population, i.e. for the 
long-time limit of Sz- Transforming the Bloch equation into the frame co- 
rotating with the driving field and computing the steady-state solution, we 
obtain 

^-('^) = 7 2 (59) 

with a linewidth 7^ = I/T2 + oj'j^2/Ti- This simple result allows spectro- 
scopic determination of all the parameters of the Bloch equation: At weak 
driving, ujr^/TiT2 <C 1, the line width is 'i/T2. This regime can be easily 
identified as the spectral line not being saturated, i.e. the height grows with 
increasing drive. In fact, the height of the resonance is Szi'yB) = lo'jiT2, 
which (knowing I/T2) allows to determine cjr. Due to the heavy filtering 
between the room-temperature driving and the cryogenic environment, this 
is not known a priory. To determine Ti , one goes to the high driving regime 
with a saturated line, i.e. a line which does not grow any more with higher 
power, u}jiy/TiT2 ^ 1 and finds a line width of iOR^jTxjT^,. With all other 
parameters known already, this allows to find T\. Using this approach is 
helpful to debug an experiment which does not work yet. Alternatively, 
real-time measurements of T\ are possible under a wide range of conditions. 
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2.2.1. How to derive the Block equation: The Block- Redfield tecknique 
We now show how to derive Bloch-hke equations from the system-bath 
models we studied before using a sequence of approximations. The Born 
approximation works if the couphng between system and bath is weak. The 
Markov approximation works if the couphng between system and bath is 
the slowest process in the system, in particular if it happens on a time scale 
longer than the correlation time of the environment. Quantitatively, we can 
put this into the motional narrowing condition 



where A is the coupling strength between the system and its environment 
and Tc the correlation time of the environment. In the case treated in eq. 
19 we would have Tc = l/wc- If this is satisfied, an averaging process over 
a time scale longer than Tc but shorter than can lead to simple evolu- 
tion equations, the so-called Bloch-Redfield equations (Argyres and Kelley, 
1964). The derivation in Ref. (Cohen- Tannoudji et al., 1992) follows this 
inspiration. We will follow the very elegant and rigorous derivation using 
projection operators as given in (Argyres and Kelley, 1964; Weiss, 1999). 
We are going to look at a quantum subsystem with an arbitrary finite 
dimensional Hilbert space, accomodating also qudit and multiple-qubit 
systems. 

As a starting point for the derivation of the Bloch-Redfield equations 
(70), one usually (Weiss, 1999) takes the Liouville equation of motion for 
the density matrix of the whole system W{t) (describing the time evolution 
of the system) 



where i^totai is the total Hamiltonian and >Ctotai the total Liouvillian of 
the whole system. This notation of the Liouvillian uses the concept of 
a superoperator. Superoperator space treats density matrices as vectors. 
Simply arrange the matrix elements in a column, and each linear operation 
on the density matrix can be written as a (super)matrix multiplication. 
Thus, the right hand side of the Liouville equation can be written as a 
single matrix products, not a commutator, where a matrix acts from the 
left and the right at the same time. Hamiltonian and Liouvillian consist of 
parts for the relevant subsystem, the reservoir and the interaction between 
these 



W{t) = [//total, M^(i)] = AotalW(t) 



(61) 




Hsys + Hj-cs + Hi 

^sys ~l~ -^res ~l~ 



(62) 
(63) 
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Hgys is the Hamiltonian which describes the quantum system (in our 
case: the qubit setup), i/res represents for the environment and Hj is the 
interaction Hamiltonian between system and bath. 

Projecting the density matrix of the whole system W{t) on the relevant 
part of the system (in our case the qubit), one finally gets the reduced 
density matrix p acting on the quantum system alone 

pit) = TvBWit) = PWit) , (64) 

so P projects out onto the quantum subsystem. As in the previous deriva- 
tion in section 1.1.1, we need to formally solve the irrelevant part of the 
Liouville equation first. Applying (1 — P), the projector on the irrelevant 
part, to eq. 61 and the obvious W = PW + (1 — P)W we get 

(1 - P)W = (1 - P)£total(l - P)W + (1 - P)£total/0. (65) 

This is an inhomogenous linear equation of motion which can be solved 
with variation of constants, yealding 

t 

{l-P)p{t) = J dt'e(i-^)^'°*-'(*"*)'(l-P)£totaip(i')+e^^~^^'^'°*^'*(l-^')W^(0). 



(66) 

Putting this result into equation (61) one gets the Nakajima-Zwanzig equa- 
tion (Nakajima, 1958; Zwanzig, 1960) 

t 

p{t) = PCtotMt)+ J dt'PCtot.ie^'-''^'-'°'^'^'-''\l - P)CtotMt') + 



+P£totaie^'"^^^'°'^'*(l - P)VF(O). (67) 

So far, all we did was fully exact. The dependence on the initial value 
of the irrelevant part of the density operator (1 — P)W{0) is dropped, if 
the projection operator is chosen appropriately - using factorizing initial 
conditions, i.e. W = (1 — P)W. A critical assesment of this assumption 
will be given in section 3.1.2. As P commutes with >Csysi one finds 

t 

p = P(£sys + Ci)p{t) + J dt'P£,e(i-^)^--' (1 - P)Cip{t'). (68) 



The reversible motion of the relevant system is described by the first (in- 
stantaneous) term of eq. (68), which contains the system Hamiltonian in 
£sys and a possible global energy shift originating from the environment 
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in RCi. The latter term can be taken into account by the redefinition 
H'g = Hs + PHj and Hj = (1 — P)Hi. The irreversibihty is given by the 
second (time-retarded) term. The integral kernel in eq. (68) still consists 
of all powers in Ci and the dynamics of the reduced density operator p of 
the relevant system depends on its own whole history. To overcome these 
difficulties in practically solving eq. (68), one has to make approximations. 
We begin by assuming that the system bath interaction is weak and restrict 
ourselves to the Born approximation, second order in Cj. This allows us 
to replace Ctotai by Cgys + Cres in the exponent. The resulting equation 
is still nonlocal in time. As it is convolutive, it can in principle be solved 
without further approximations (Loss and DiVincenzo, 2003). To proceed 
to the more convenient Bloch-Redfield limit, we remove the memory firstly 
by propagating p{t') forward to p{t). In principle, this would require solving 
the whole equation first and not be helpful. In our case, however, we can 
observe that the other term in the integral — the kernel of the equation — is 
essentially a bath correlation function which only contributes at t — t' < Tc- 
Using the motional narrowing condition eq. 60, we see that the system 
is unlikely to interact with the environment in that period and we can 
replace the evolution of p with the free eveloution, p{t') = e^'^^'sC*-* ^ p{t). 
After this step, the equation is local in time, but the coefficients are still 
time-dependent. Now we flip the integration variable t' ^ t — t' and then 
use the motional narrowing condition again to send the upper limit of the 
integral to infinity, realizing that at such large time differences the kernel 
will hardly contribute anyway. We end up with the Bloch-Redfield equation 

oo 

p{t) = P{Csy, + Ci)p{t) + J dt'P£,e(i-^)(^-+^-)*'(l - P)Cip{t). (69) 



The Bloch-Redfield equation is of Markovian form, however, by properly 
using the free time evolution of the system (back-propagation), they take 
into account all bath correlations which are relevant within the Born ap- 
proximation (Hartmann et al., 2000). In (Hartmann et al., 2000), it has 
also been shown that in the bosonic case the Bloch-Redfield theory is 
numerically equivalent to the path-integral method. 

The resulting Bloch-Redfield equations for the reduced density matrix 
p in the eigenstate basis of Hgys then read (Weiss, 1999) 

Pnm{t) = —i^nmPnm{t) — RnmklPkljt) i (70) 

where Rnmki ^ire the elements of the Redfield tensor and the pnm are the 
elements of the reduced density matrix. 
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The Redfield tensor has the form (Weiss, 1999; Blum, 1996) 

Rnmke = ^em X! '^nrlk + ^nk X! '^irrm ~ ^tmnk ~ ^emnk' i'^^) 

r r 

The rates entering the Redfield tensor elements are given by the follow- 
ing Golden-Rule expressions (Weiss, 1999; Blum, 1996) 

CXD 

rJifc = dt e-^'^-^\Hi,Ut)Hi,nkm (72) 



oo 

rLL = ?^"'/ dt e-'^'-'{HjMmi,nk{t)) , (73) 



where Hj appears in the interaction representation 

Hj{t) = eM^Hrest/h) Hj exp{-iH,^,t/h). (74) 

Unk is defined as Unk = {En — Ek)/h. In a two-state system, the coefficients 
i, m, n and k stand for either -|- or — representing the upper and lower 

eigenstates. The possible values of ujnk in a TSS are = uj = 0, 

W-i = ^ and u) |_ = — where E is the energy splitting between the 

two charge eigenstates with E = Ve^ + A^. Now we apply the secular 
approximation, which again refers to weak damping, to discard many rates 
in the Redfield tensor as irrelevant. The details of this approximation are 
most transparent in the multi-level case and will be discussed in more detail 
in section 4.0.4. In the TSS case, the secular approximation holds whenever 
the Born approximation holds. After the secular approximation, the Bloch- 
Redfield equation coincides with the Bloch equation with 

1/^1 = Rnnnn = + R = r-++- -|- r_++_ (75) 

n 

I/T2 = Re{Rnmnm) = Re(i?+_+_) = Re(i?_+_+) 

= Re(r+„„+ + r_++_ + r + r++++ - r__++ - r++„„). 

= ^ + ;^ (76) 

Here, we have introduced the dephasing rate T^^ . The relaxation rate is 
given by the time evolution of the diagonal elements, and the dephasing 
rate by the off- diagonal elements of the reduced density matrix p. 

The factor of two in the formula connecting I/T2 and 1/Ti appears to 
be counterintuitive, as we would expect that energy relaxation definitely 
also leads to dephasing, without additional factors. This physical picture is 



SUPERCONDUCTING QUBITS II: DECOHERENCE 23 



S.e, = '.^^^-^^ = tanh ( ^ ) (77) 



also correct, but one has to take into account that there are two channels for 
dephasing — clockwise and counterclockwise precession — which need to be 
added. In fact, this is the reason why the same factor of two appears in the 
positivity condition for the density matrix, see section 2.1.3. Another view 
is to interpret the diagonal matrix elements as classical probabilities, the 
absolute square of a eigenfunctions of the Hamiltionian, iV'iP; whereas the 
off-diagonal terms constitute amplitudes, V'2^i- Being squares, probabilities 
decay twice as fast as amplitudes. This point will be discussed further later 
on in the context of multi-level decoherence, eq. 106. 

The imaginary part of the Redfield tensor elements that are relevant for 

the dephasing rate ^(i?-! i ) provides a renormalization of the coherent 

oscillation frequency a; _| , 5uj^ = 9(r^ h + T hH )• If the renormaliza- 
tion of the oscillation frequency gets larger than the oscillation frequency 
itself, the Bloch-Redfield approach with its weak-coupling approximations 
does not work anymore. By this, we have a direct criterion for the validity 
of the calculation. 

Finally, the stationary population is given by 

= tanh — — - 

+ r+__+ \2kBT, 

where in the last step we have used the property of detailed balance 

P — "T „-l^mn/kBT /rjo) 

nmmn — mnnm'^ \ ' "/ 

which holds for any heat bath in thermal equilibrium and is derived e.g. in 
References (Weiss, 1999; Ingold, 1998; Callen and Welton, 1951). 

A different kind of derivation with the help of Keldysh diagrams for 
the specific case of an single-electron transistor (SET) can be found in the 
Appendix of Ref. (Makhlin et al., 2001). 

Very recent results (Gutmann and Wilhelm, 2006; Thorwart et al., 2005) 
confirm that without the secular approximation, Bloch-Redfield theory pre- 
serves complete positivity only in the pure dephasing case (with vanishing 
coupling A = between the qubit states). In all other cases, complete 
positivity is violated at short time scales. Thus only in the pure dephas- 
ing regime is the Markovian master equation of Lindblad form (Lindblad, 
1976) as typically postulated in mathematical physics. In all other cases the 
Lindblad theorem does not apply. This is not an argument against Bloch- 
Redfield — the Markovian shape has been obtained as an approximation 
which coarse-grains time, i.e. it is not supposed to be valid on short time 
intervals. Rather one has to question the generality of the Markov approx- 
imation (Lidar et al., 2004) at low temperature. Note, that in some cases 
the violation of positivity persists and one has to resort to more elaborate 
tools for consistent results (Thorwart et al., 2005) 
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2.2.2. Rates for the Spin-Boson model and their physical meaning 
This technique is readily apphed to the spin boson Hamiltonian eq. (32). 
The structure of the golden rule rates eqs. (72 and 73) become rather 
transparent — the matrix elements of the interaction taken in the energy 
eigenbasis measure symmetries and selection rules whereas the time integral 
essentially leads to energy conservation. 

In particular, we can identify the energy relaxation rate 

^ = ^S(E). (79) 

The interpretation of this rate is straightforward — the system has to make 
a transition, exchanging energy E with the environment using a single 
Boson. The factor S{E) = J[E){n[E) + 1 + n[E)) captures the density of 
Boson states J{E) and the sum of the rates for emission proportional to 
n{E) + 1 and absorption proportional to n{E) of a Boson. Here, n{E) is 
the Bose function. The prefactor is the squared cosine of the angle between 
the coupling to the noise and the qubit Hamiltonian, i.e. it is maximum 
if — in the basis of qubit eigenstates — the bath couples to the qubit in 
a fully off-diagonal way. This is reminiscent of the standard square of the 
transition matrix element in Fermi's golden rule. 
The flip-less contribution to T2 reads 

(80) 

It accounts for the dephasing processes which do not involve a transition of 
the qubit. Hence, they exchange zero energy with the environment and 5(0) 
enters. The prefactor measures which fraction of the total environmental 
noise leads to fluctuations of the energy splitting, i.e., it is complemetary 
to the transition matrix element in Ti — the component of the noise 
diagonal in the basis of energy eigenstates leads to pure dephasing. The 
zero frequency argument is a consequence of the Markov approximation. 
More physically, it can be understood as a limiting procedure involving the 
duration of the experiment, which converges to S{Q) under the motional 
narrowing condition. Details of this procedure and its limitations will be 
discussed in the next section. 
Finally, the energy shift 

SE = ^,vj d.^^, (81) 

where V denotes the Cauchy mean value, is analogous to the energy shift 
in second order perturbation theory, which collects all processes in which a 
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virtual Boson is emitted and reabsorbed, i.e. no trace is left in the environ- 
ment. Again, the prefactor ensures that the qubit makes a virtual transition 
during these processes. For the Ohmic case, we find 



provided that ^ E. Thus, the energy shift explicitly depends on the 
ultraviolet cutoff. In fact, 6E ~ E would be an indicator for the breakdown 
of the Born approximation. Thus, we can identify two criteria for the valid- 
ity of this approximation, a <^ 1 and a\og{uJc/E) <^ 1. The latter is more 
confining, i.e. even if the first one is satisfied, the latter one can be violated. 
Note that in some parts of the open quantum systems literature, the justi- 
fication and introduction of this ultraviolet cutoff is discussed extensively. 
The spectral densities we have computed so far in the previous sections have 
always had an intrinsic ultraviolett cutoff, e.g. the pure reactive response 
of electromagnetic circuits at high frequencies. 

2.3. ENGINEERING DECOHERENCE 

The picture of decoherence we have at the moment apparently allows to 
engineer the decoherence properties — which we initially percieved as some- 
thing deep and fundamental — using a limited set of formulae, eqs. 79, 80 
and 42, see Refs. (van der Wal et al., 2003; Makhhn et al., 2001) these 
equations have been applied to designing the circuitry around quantum 
bits. This is, however, not the end of the story. After this process had been 
mastered to sufficient degree, decoherence turned out to be limited by more 
intrinsic phenomena, and by phenomena not satisfactorily described by the 
Bloch-Redfield technique. This will be the topic of the next section. 

3. Beyond Bloch-Redfield 

It is quite surprising that a theory such as Bloch-Redfield, which contains 
a Markov approximation, works so well at the low temperatures at which 
superconducting qubits are operated, even though correlation functions at 
low temperatures decay very slowly and can have significant power-law 
time tails. The main reason for this is the motional narrowing condition 
mentioned above, which essentially states that a very severe Born approx- 
imation, making the system-bath interaction the lowest energy /longest 
time in the system, will also satisfy that condition. This is analogous to 
the textbook derivation of Fermi's golden rule (Cohen- Tannoudji et al., 
1992; Sakurai, 1967), where the perturbative interaction is supposed to be 




(82) 
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the slowest process involved. In this section, we are going to outline the 
limitations of this approach by comparing to practical alternatives. 

Before proceeding we would also like to briefly comment on the general 
problem of characterizing the environment in an open quantum system. The 
most general environment is usually assumed to induce a completely posi- 
tive linear map (or "quantum operation") on the reduced density matrix. 
The most general form of such a map is known as the Krauss operator- 
sum representation, although such a representation is not unique, even 
for a given microscopic system-bath model like the one considered here. 
A continuous-time master equation equivalent to a given Krauss map is 
provided by the Lindblad equation, but the form of the Lindblad equa- 
tion is again not unique. The Lindblad equation gives the most general 
form of an equation of motion for the reduced density matrix that assures 
complete positivity and conserves the trace; however, the Marvok and Born 
approximations are often needed to construct the specific Lindblad equation 
corresponding to a given microscopic model. The Markov approximation 
is a further additional simplification, rendering the dynamics to that of 
a semigroup. A semigroup lacks an inverse, in accordance with the un- 
derlying time-irreversibility of an open system. However, like the unitary 
group dynamics of a closed system, the semigroup elements can be gener- 
ated by exponentials of non-Hermitian " Hamiltonians" , greatly simplifying 
the analysis. The Bloch-Redfield master equation also has a form similar 
to that of the Lindblad equation, but there is one important difference: 
Bloch-Redfield equation does not satisfy complete positivity for all values 
of the diagonal and off-diagonal relaxation parameters. If these parameters 
are calculated microscopically (or are obtained empirically), then complete 
positivity will automatically be satisfied, and the Bloch-Redfield equation 
will be equivalent to the Lindblad equation. Otherwise inequalities have to 
be satisfied by the parameters in order to guarantee complete positivity. 

3.1. PURE DEPHASING AND THE INDEPENDENT BOSON MODEL 

We start from the special case A = of the spin-Boson model, also known 
as the independent Boson model (Mahan, 2000). We will discuss, how 
this special case can be solved exactly for a variety of initial conditions. 
Restricting the analysis to this case is a loss of generality. In particular, as 
the qubit part of the Hamiltonian commutes with the system-bath coupling, 
it cannot induce transitions between the qubit eigenstates. Thus 1/Ti = 
to all orders as confirmed by eq. 79 and I/T2 = S{0) following eq. 80. Still, 
it allows to gain insight into a number of phenomena and the validy of the 
standard approximations. Moreover, the results of this section have been 
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confirmed based on a perturbative diagonalization scheme valid for gap or 
super-olimic environmental spectra (Wilhelm, 2003). 

3.1.1. Exact propagator 

As the qubit and the qubit-bath coupling commute, we can construct the 
exact propagator of the system. We go into the interaction picture. The 
system-bath coupling Hamiltonian then reads 

HsBit) = Aj(a.e-*'^^* + a\e'^^'). (83) 

i 

The commutator of this Hamiltonian with itself taken at a different time is 
a c-number. Consequently, up to an irrelevant global phase, we can drop the 
time-ordering operator T in the propagator (Sakurai, 1967; Mahan, 2000) 
and find 

U(t,t') = Texp(-^l'dt'HsB{t')^ (84) 

= -p E ^ (4 (»■"•"-'■'-)-".('-"•"-'''-'))) ■ 

In order to work with this propagator, it is helpful to reexpress it using 
shift operators -Dj(ai) = exp(aa''' — a* a) as 



e^"^ - l) . (85) 



This propagator can be readily used to compute observables. The main tech- 
nical step remains to trace over the bath using an appropriate initial state. 
The standard choice, also used for the derivation of the Bloch-Redfield equa- 
tion, is the factorized initial condition with the bath in thermal equilbrium, 
i.e. the initial density matrix 

where we use the partition function Z (Landau and Lifshitz, 1984). The 
expectation value of the displacement operator between number states is 
(n| D{a) \n) = e"^^"'"^^)!'^! We start in an arbitrary pure initial state of 
the qubit 

p, = |V)(V|, IV) =cos^|0) + sin^e'<^|l). (87) 

Using these two expressions, we can compute the exact reduced density 
matrix, expressed through the three spin projections 

{a^){t) = sin6'cos(£;t + (/))e--^/W 
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{ay){t) = sm9sm{Et + (t>)e-^f^^^ (89) 
{az){t) = cose (90) 

where we have introduced the exponent of the envelope for factorized initial 
conditions, 

Kf{t) = J ^S{lo){1- cos ujt) (91) 

which coincides with the second temporal integral of the semiclassical cor- 
relation function S{t), see eq. 18. What does this expression show to us? 
At short times, we always have Kf{t) oc y / dujS{uj)^ which is an integral 
dominated by large frequencies and thus usually depends on the cutoff of 
S{u). At long times, it is instructive to rewrite this as 



Kf{t) = t J duj6^{t)S{uj) (92) 

where we have introduced 6co{t) = 2 ^^'^Ji^^^ , which approaches 6{lo) as 
t — > oo. Performing this limit more carefully, we can do an asymptotic 
long-time expansion. Long refers to the internal time scales of the noise, 
i.e. the reciprocal of the internal frequency scales of S{uj), including h/kT, 
uj~^. The expansion reads 

Kf{t) = -t/T2 + log VF + 0{l/t) (93) 

with I/T2 = 5(0) as in the Bloch-Redfield result and logvp = Vj^S{u;). 
Here, V is the Cauchy mean value regularizing the singularity at a; = 0. To 
highlight the meaning ol vp, the visibility for factorized initial conditions, 
we plug this expansion into eq. 88 and see that {ax) (t) = vp sm9cos{Et + 
(/))e~*/"^^''~'^^^/*^. Thus, a long-time observer of the full dynamics sees expo- 
nential decay on a time scale T2 which coincides with the Bloch-Redfield 
result for the pure dephasing situation, but with an overall reduction of 
amplitude by a factor v < 1. This is an intrisic loss of visibility (Vion et al., 
2002; Simmonds et al., 2004). Several experiments have reported a loss of 
visibility, to which this may be a contribution. Note that by improving 
detection schemes, several other sources of reduced visibility have been 
eliminated (Lupascu et al., 2004; Wallraff et al., 2005). 

This result allows a critical assessment of the Born-Markov approxi- 
mation we used in the derivation of the Bloch-Redfield equation. It fails 
to predict the short-time dynamics — which was to be expected as the 
Markov approximation is essentially a long-time limit. In the long time 
limit, the exponential shape of the decay envelope and its time constant are 
predicted correctly, there are no higher-order corrections to T2 at the pure 
dephasing point. The value of T2 changes at finite A, see Refs. (Leggett 
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et al, 1987; Weiss, 1999; Wilhelm, 2003). A further description of those 
results would however be far beyond the scope of this chapter and can 
be found in ref (Wilhelm, 2003). Finally, we can see how both short and 
long-time dynamics are related: the short-time (non-Markovian) dynamics 
leaves a trace in the long-time limit, namely a drop of visibility. 

We now give examples for this result. In the Ohmic case J{co) = aLoe~'^l^'' 
at T = 0. Hence, we can right away compute Kf{t) and obtain Kf{t) = 
^ log(l -|- (wci)^) by a single time integral. In agreement with the formula 
for T2, see eqs. 76, 80, the resulting decay does not have an exponential 
component at long time but keeps decaying as a power law, indicating 
vanishing visibility. 

At finite temperature, the computation follows the same idea but leads 
to a more complicated result. We give the expression from Ref. (Gorlich 
and Weiss, 1988) for a general power-law bath Jq{uj) = aqijj'^uj].~'^e~^/^'' ^ 

Kf{t) = 2Re |a,r(g _ 1) j^i _ (1 + iu;,t)^-^ + (^^^ x (94) 
r w ^^ ^ ikTt\ ( ^ ikTkT\]\] 

X 2C{s-i,n)-c^s-i,n + ^j -C^q-i,n-^—j^jj 

where we have introduced ^1 = 1 + kBT/hujQ and the generalized Riemann 
zeta function, see (Abramowittz and Stegun, 1965) for the definition and the 
mathematical properties used in this subsection. This exact result allows to 
analyze and quantify the decay envelope by computing the main parameters 
of the decay, vp and I/T2. We will restrict ourselves to the scaling limit, 
iOc ^ l/t,kT. For the Ohmic case, q = 1, we obtain at finite temperature 
I/T2 = 2akT/h and vp = {kT /ojc)"'. This result is readily understood. The 
form of T2 accounts for the fact that an Ohmic model has low-frequency 
noise which is purely thermal in nature. The visibility drops with growing 
u)c indicating that if we keep adding high frequency modes they all con- 
tribute to lost visibility. It is less intuitive that vp drops with lowering the 
temperature, as lowering the temperature generally reduces the noise. This 
has to be discussed together with the l/T2-term, remembering that I/T2 is 
the leading and vp only the sub- leading order of the long time expansing 
eq. 93: At very low temperatures, the crossover to the exponential long- 
time decay starts later and the contribution of non-exponential short time 
dynamics gains in relative significance. Indeed, at any given time, the total 
amplitude gets enhanced by lowering the temperature. 

In order to emphasize these general observations, let us investigate the 
super-Ohmic case with g > 3. Such spectral functions can be realized in 
electronic circuits by RC-series shunts (Robertson et al., 2005), they also 
play a significant role in describing phonons. For (7 > 3, the exponential 
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component vanishes, I/T2 = and Vg = exp[—2aqT{q — l)]. Thus, we obtain 
a massive loss of visibiHty but no exponential envelope at all. This highlights 
the fact that v and I/T2 are to be considered independent quantifiers of 
non-Markovian decoherence and that the latter accounts for environmental 
modes of relatively low frequency whereas v is mostly influenced by the fast 
modes between the qubit frequency and the cutoff. 

Before outlining an actual microscopic scenario, we generalize the initial 
conditions of our calculation. 

3.1.2. Decoherence for non-factorzing initial conditions 

Our propagator, eq. 85, is exact and can be applied to any initial density 

matrix. We start from an initial wave function 



IV') = |0) l[^D{z^/2X,huj,) |0), + |1) l[^D{zj/2X,hco,) |0), (95) 



where we have introduced sets of dimensionless coefficients z^ .It would 
be straightforward to introduce 6 and (j), which we will stay away from 
in order to keep the notation transparent. The factorized initial condition 
corrsponds to z^^^ = 0. 

This structure has been chosen in order to be able to obtain analytical 
results, using the structure of the propagator expressed in displacement 
operators, eq. 85 and the multiplication rules for these operators (Walls 
and Milburn, 1994). Note that the choice of coherent states to entangle the 
qubit with is not a severe restriction. It has been shown in quantum optics in 
phase space, that essentially each density matrix of an harmonic oscillator 
can be decomposed into coherent states using the Wigner or Glauber P 
phase space representations, see e.g. (Schleich, 2001). Physically, the initial 
state eq. 95 corresponds to the qubit being in a superposition of two dressed 
states. Of specific significance is the initial condition which minimizes the 
sytem bath-interaction in the Hamiltonian eq. 32, nameley zf = —z~^ = 
-1. 

We can again compute all three spin projections of the qubit. The 
essence of the decoherence behavior is captured in the symmetric initial 
state, z'- = —zj for all i 




(96) 



very similar to eq. 88 in the factorized case, but now with 



K{t) 



1 dii} V 

- J{^u:)[{u{u:) + lf+v\u:) + l 

2 Jo LO^ L 



2(1 + u{uj) cosujt + v{uj) shiLot)]) 
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where we have taken a continuum hmit replacing the complex numbers 
by the real function u{lo) + iv{Lv). This form connects to the factorized case 
by setting u = v = 0. For any other choice of u and v, the initial conditions 
are entangled. 

We can make a few basic observations using this formula: The initial 
amplitude e~^^^^ is controlled through 



thus for any initial condition which is more than marginally entangled 
(meaning that the integral is nonzero), the initial amplitude is smaller than 
unity. On the other hand, the time-dependence of K{t) can be completely 
eliminated by chosing an initial condition u = —l,v = 0. This condition 
minimizes the system-bath part of the total energy in the sense of variation 
with respect to u and v. This choice of initial state also minimizes the total 
energy if the oscillators are predominantly at high frequency, whereas for 
the global minimum one would rather chose a factorized state for the low- 
frequency oscillators. Physically, this corresponds to an optimally dressed 
state of the qubit surrounded by an oscillator dressing cloud. The overlap 
of these clouds reduces the amplitude from the very beginning but stays 
constant, such that the long-time visibility 



is maximum. Note that this reduces to the result for vp for u = v = 0. 

What can we learn from these results? We appreciate that initial condi- 
tions have a significant and observable effect on the decoherence of a single 
qubit. The choice of the physically appropriate initial condition is rather 
subtle and depends on the experiment and environment under considera- 
tion. A free induction decay experiment as described here does usually not 
start out of the blue. It is launched using a sequence of preparation pulses 
taking the state from a low temperature thermal equilibrium to the desired 
initial polarization of the qubit. Thus, from an initial equilbrium state (for 
some convenient setting of the qubit Hamiltonian) , the fast preparation 
sequence initiates nonequilibrium correlations thus shaping u and v. Fur- 
thermore, if the interaction to the environment is tunable such as in the 
case of the detectors discussed previously in section 2.1.2, the initial condi- 
tion interpolates between factorized (rapid switching of the qubit-detector 
coupling) and equilibrium (adiabatic switching). 

At this point, we can draw conclusions about the microscopic mechanism 
of the loss of visibility and other short-time decoherence dynamics. The 
picture is rooted on the observation that the ground state of the coupled 




(97) 




(98) 



32 F.K. WILHELM, M.J. STORCZ, U. HARTMANN, AND M. GELLER 



system is a dressed state. On the one hand, as described above, the overlap 
of the dressing clouds reduces the final visibility. On the other hand, for 
nonequilibrium initial conditions such as the factorized one, there is extra 
energy stored in the system compared to the dressed ground state. This 
energy gets redistributed while the dressing cloud is forming, making it 
possible for an excitation in the environment with an extra energy 5E to 
be created leading to a virtual intermediate state, followed by another exci- 
tation relaxing, thus releasing the energy 6E again. It is crucial that this is 
another excitation as only processes which leave a trace in the environment 
lead to qubit dephasing. Higher-order processes creating and relaxing the 
same virtual excitation only lead to renormalization effects such as the 
Lamb shift, see eq. 81. This explains why the loss of visibility is minimal 
for dressed initial conditions, where no surplus excitations are present. 

The Bloch-Redfield technique is a simple and versatile tool which makes 
good predictions of decoherence rates at low damping. At higher damp- 
ing, these rates are mostly joined by renormalization effects extending the 
Lamb shift in eq. 81, see Refs. (Leggett et al, 1987; Weiss, 1999; Wil- 
helm, 2003). However, there is more to decoherence than a rate for accu- 
rate predictions of coherence amplitudes as a function of time, one has 
to take the non-exponential effects into account and go beyond Bloch- 
Redfield. Other approaches can be applied to this system such as rigor- 
ous (Born but not Markov) perturbation theory (Loss and DiVincenzo, 
2003), path- integral techniques (Leggett et al., 1987), (Weiss, 1999), and 
renormalization schemes (Kehrein and Mielke, 1998). 

Note, that these conclusions all address free induction decay. There is 
little indication on the quality of the Bloch-Redfield theory in the presence 
of pulsed driving. 

3.1.3. 1/f noise 

In the previous sections we have explored options how to engineer deco- 
herence by influencing the spectral function J{io) e.g. working with the 
electromagnetic environment. This has helped to optimize supercondcuting 
qubit setups to a great deal, down to the level where the noise intrinsic to 
the material plays a role. In superconductors, electronic excitations are 
gapped (Tinkham, 1996) and the electron phonon interaction is weak due 
to the inversion symmetry of the underlying crystal everywhere except 
close to the junctions (loffe et al., 2004). The most prominent source of 
intrinsic decoherence is thus 1/f noise. 1/f noise - noise whose spectral 
function behaves following S{u)) oc 1/uj, is ubiquitous in solid-state systems. 
This spectrum is very special as all the integrals in our discussion up to 
now would diverge for that spectrum. 1/f typically occurs due to slowly 
moving defects in strongly disordered materials. In Josephson devices, there 
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is strong evidence for 1// noise of gate charge, magnetic flux, and critical 
current, leading to a variety of noise coupling operators (see Ref. (Harlingen 
et al., 1988) for an overview). Even though there does not appear to be a 
fully universal origin, a "standard" model of 1// noise has been identified 
(Dutta and Horn, 1981; Weissman, 1988): The fundamental unit are two- 
state fluctuators, i.e. two state systems which couple to the device under 
consideration and which couple to an external heat bath making them 
jump between two positions. The switching process consists of uncorrelated 
switching events, i.e. the distribution of times between these switches is 
Poissoinan. If we label the mean time between switches as r, the spectral 
function of this process is Sutn = Sq i^l-I^-2 ■ This phenomenon alone is 
called random telegraph noise (RTN). Superimposing such fluctuators with 
a flat distribution of switching times leads to a total noise spectrum pro- 
portional to 1//. Nevertheless, the model stays different from an oscillator 
bath. The underlying thermodynamic limit is usually not reached as it is 
approached more slowly: Even a few fluctuators resemble 1/f noise within 
the accuracy of a direct noise measurement. Moreover, as we are interested 
in very small devices such as qubits, only a few fluctuators are effective and 
experiments can often resolve them directly (Wakai and van Harlingen, 
1987). Another way to see this is to realize that the RTN spectrum is 
highly non-Gaussian: A two - state distribution can simply not be fitted 
by a single Gaussian, all its higher cumulants of distribution are relevant. 
This non-Gaussian component only vanishes slowly when we increase the 
system size and is significant for the case of qubits. 

A number of studies of models taking this aspect into account have 
been published (Paladino et al., 2002; Grishin et al., 2005; Shnirman et al., 
2005; Galperin et al., 2003; Faoro et al., 2005; de Sousa et al., 2005). A 
highly simplified version is to still take the Gaussian assumption but realize 
that there is always a slowest fluctuator, thus the integrals in Kf{t) can be 
cut off at some frequency ujir at the infrared (low frequency) end of the 
spectrum, i.e. using the spectral function 



with 9 the Heaviside unit step function, we approximately find (Cottet, 
2002; Martinis et al., 2003; Shnirman et al., 2002) 




•2 



(99) 



(100) 



so we find the Gaussian decay typical for short times - short on the scale 
of the correlation time of the environment, which is long as the spectrum 
is dominated by low frequencies - with a logarithmic correction. 
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At the moment, forefront research works at understanding more detailed 
models of 1// noise and understand the connection between the strong 
dephasing and a possible related relaxation mechanism at high frequencies. 
On the other hand, experiments work with materials to avoid l//-noise at 
its source. Generally, slow noise up to a certain level can be tolerated using 
refocusing techniques such as simple echo or the Carr-Purcell-Gill-Meiboom 
pulse sequence (Gutmann et al., 2005; Carr and Purcell, 1954; Faoro and 
Viola, 2004; Falci et al., 2004; Shiokawa and Lidar, 2004; Bertet et al., 
2005b), the power and potential of which has been demonstrated both 
experimentally and theoretically. 

4. Decoherence in coupled qubits 

To conclude, we want to outline how to go beyond a single to multiple qubits 
and identifty the underlying challenges. On that level, much less is known 
both theoretically and experimentally. The variety of physically relevant 
Hamiltonians is larger. One extreme case is fully uncorrelated noise, e.g. 
originating from effects in the junctions or qubit-specific Hamiltonians, 

H = Hqi + Hq2 + Hqq + Hqibi + Hbi + Hq2B2 + Hb2 (Id) 

this is simply the sum of two single-qubit decoherence Hamiltonians in 
distinct Hilbert spaces, consisting of qubit Hamiltonians -ffgi, i = 1,2, 
baths Hbi, qubit-bath interaction HgiBi all interacting via a qubit-qubit 
interaction Hqq alone. The other extreme case is collective noise, e.g. 
long-wavelength ambient fluctuations or noise shared control lines. This 
is described by 

H = Hqi + Hq2 + Hqq + Hqib + Hq2B + Hb (102) 

where both qubits talk to a single bath. The distinction of baths may seem 
artifical, as this is a special case of Hamiltonian 101: What we really mean 
is that in the interaction picture there is a significant correlation between 
baths {HQiB{t)HQ2B{t')) / 0. Note, that intermediate cases between these, 
a partially correlated model (Storcz et al., 2005a), can be identified in the 
context of quantum dots. 

4.0.4. The secular approximation 

What does it take to study decoherence here or in other multilevel systems? 
Basically we can follow all the steps through the derivation of the Bloch- 
Redfield equation given in section 2.2.1 up to eq. 74 until we solve the 
equation. There, we have already mentioned the secular approximation 
without explaining its details. 
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The essence of the secular approximation is the separation of time scales. 
We go back to the interaction representation of equation 70, leading to 

= E ^n^fc^e^^--— )VL. (103) 

kl 

As the Bloch-Redfield equation is based on a Born approximation, we can 
expect \Rijki\ ^ ^mn for all coefficients i,j,k,l,m,n with m ^ n. 

In the secular limit, this also holds true for most frequency splittings 

|w„m -UJkl\ > \Rnmkl\ (104) 

besides the inevitable exceptions n = m, k = I, and n = k and m = I. 
Whenever condition eq. 104 is satisfied, the time evolution induced by Rnmki 
is certainly slower than the precession with uJnm — ^ki and averages out 
quickly, hence it can be dropped. So the only remaining rates are the cases 
just mentioned: 

For n = k,m = I, we have to keep Rnmnm- This rate is the dephasing 
rate for the transition between levels nm, see eq. 80. These rates depend 
on the pair of levels we chose and in general they will all be different for 
different choices of n and m, leading to N{N — l)/2 different T2-rates for 
an A^-level system. 

For n = m and k = I. The set of these terms splits off from the rest of 
the equation, i.e. the diagonal terms of the density matrix (in the eigenstate 
basis) decay independent from the off-diagonal terms and obey the following 
set of equations 

Pn — ^^(-fm-Tm— >}2 n^m) (105) 

n 

which is analogous to the Pauli master equation for classical probabilities. 
We have identified P„ = pnm the classical probability and the transition 
rates r„^m = Rnnmm- Equation 79 can be solved by Laplace transform, 
where it reduced to a matrix inversion. This leads to different inde- 
pendent energy relaxation channels whose rates are the eigenvalues of the 
matrix form of the right hand side of eq. 105. One of these eigenvalues is 
always zero representing stable thermodynamic equilbrium which does not 
decay. In the two-state case, this leads us to one nonvanishing Ti rate rep- 
resenting the only nonzero eigenvalue, given by eq. 79. The rates generally 
obey the positivity constraint 

^ ^ Rnnnn ^ 2 ^ ^ Rnmnm (106) 
n n^tm 

the left hand side being the trace of the relaxation matrix, i.e. the sum over 
all Ti-type rates and the ricght hand side being the sum over dephasing 
rates. This reduces to T2 < 2Ti in a two-state system. 
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In the opposite case, the case of a approximate Liouvillian degeneracy, 
we find a pair of frequencies \u;nm — ^ki\ ^ \Rnmki\ which do not obey 
the conditions mentioned in the previous paragraph, such that the secu- 
lar approximation does not apply to this set of levels and Rijki must be 
kept. In that case, the Bloch-Redfield equation can still be diagonalized 
numerically, identifying the relvant modes of decay (van Kampen, 1997). 
Note, that Liouvillian degeneracies can appear in non-degenerate systems, 
promintently the single harmonic oscillator. One practical example for this 
issue is intermediate-temperature cavity QED (Rau et al., 2004) 

These concepts already found some application in the theoretical liter- 
ature. We just mention the main results here. After the pioneering work 
(Governale et al., 2001), it was realized that the high number of rates makes 
the results difficult to analyze and the performance of quantum gates should 
be analyzed directly (Thorwart and Hanggi, 2002; Storcz and Wilhelm, 
2003; Wilhelm et al., 2003). A key result is that (only) the correlated noise 
model, eq. 102 permits to use symmetries and encoding into decoherence 
free subspaces to protect coherence (Storcz and Wilhelm, 2003; Wilhelm 
et al., 2003; Storcz et al., 2005c), where deviations from perfect symmetry 
are of relatively low impact (Storcz et al., 2005b). 

5. Summary 

In summary, we have provided an introduction to standard methods in 
decoherence theory as they are applied to superconducting qubits. Many of 
the tools and results are more general and can be applied to other damped 
two-state-systems as well. We see that parts of the theory of decoherence — 
in particular the part on electromagnetic environments and Bloch-Redfield- 
Theory — are really well established by now, only opening the view on 
more subtle problems connected to memory effects and the interplay of 
decoherence and control. 
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